On optimal values of a for the analytic Hartree-Fock-Slater method 
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We have examined the performance of the analytic Hartree-Fock-Slater (HFS) method for various 
a values and empiricaly determined the optimal a value by minimizing the mean absolute error 
(MAE) in atomization energies of the G2 set of molecules. At the optimal a the HFS method's 
performance is far superior with the MAE of 14 kcal/mol than that of the local density approximation 
(MAE ~ 36 kcal/mol) or the Hartree-Fock theory (MAE ~ 78 kcal/mol). The HFS exchange 
functional with a — 0.7091 performs significantly better than the Kohn-Sham exchange functional 
for equally weighted atoms H-Kr. We speculate that use of this single a value may be useful in 
parametrization of empirical exchange-correlation functionals. 
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The problems with numerical integration in quantum 
mechanical calculations are well knowni^ Thus, all ab 
initio electronic structure calculations can be divided 
into two classes. In the first class, historically called 
ab initio, all quantum mechanical matrix elements are 
computed to machine precision^ The second class of 
electronic structure calculations require numerical inte- 
gration and machine-precision matrix elements are to- 
tally impractical, and thus, except for atoms, machine- 
precision energies are out of the question. Until very 
recently this second class included almost all density- 
functional calculations. 4 Recently, fully analytic Hartree- 
Fock-Slater— (HFS) variant of the density functional the- 
ory (DFT) was implemented using Gaussian basis setsA 
This approach employs fitting of the potential to inte- 
grable functional form, rather than by fitting or integra- 
tion on numerical grid. The technique is computationally 
very efficient in comparison with the grid-based imple- 
mentation and provides smooth potential energy surfaces 
and exact energy gradients^ We have recently extended 
this scheme to allow for the atom-dependent exchange 
parameters a that scale the exchange potential by means 
of a muffin-tin (MT)-like approach^ In our method ma- 
trix elements arc computed to machine accuracy. Fur- 
ther, in contrast to earlier MT implementation, here the 
energy is both meaningful and stationary. One can re- 
quire that atoms dissociate into their exact experimental 
rather than approximate Hartree-Fock (HF) electronic 
energies. This approach^ when applied to a standard set 
of molecules that are used in performance tests of DFT 
models yields results that are intermediate between ei- 
ther the local density approximation (LDA) or the HF 
appoxumation and more sophisticated hybrid or general- 
ized gradient approximations (GGA)A 



In the HFSS model, the nonlocal exchange potential in 
the Hartree-Fock method is replaced by a local exchange 



potential that is given by 

*w=-§(f) 1/ V /3 - a) 

Here, the parameter a, called Slater's statistical exchange 
parameter, is unity. Similar expression for the exchange 
energy of the homogeneous electron gas was obtained ear- 
lier by DiraciS Later, Gaspar— and Kohn-Sham-A (GKS) 
obtained the value of 2/3 for a by variationally minimiz- 
ing the total energy functional. In the following years, 
a was taken purely as an adjustable parameter to ob- 
tain desired atomic propertiesi 12 i 13 i 14 i 15 The firsti^ HFS 
calculations with meaningful numerically integrated to- 
tal energies used a uniform a value of 0.7. Since then 
the HFS method has come to mean this a value. Later, 
the electronic structure calculations using the LDA by 
showed that the LDA give similar— but not superior- 
binding energies to the HFS method. Several studies 
since then have shown that the LDA has a general ten- 
dency to overbindfi^ 

HF theory being analytic allows cheap geometry opti- 
mization despite its N cost. In an analytic method one 
optimizes tens of linear-combination-of-atomic-orbital 
parameters per atom, rather than hundreds of plane- 
waves per pseudoatom, or thousands of numerical inte- 
gration points per all-electron atom. With or without the 
MT-likc advance, an N 3 analytic method might prove 
to be a practical geometry-optimization tool if appropri- 
ate choice(s) of the exchange parameter(s) is(are) made. 
In this article we asses the performance of analytic HF 
model for the GKS and the Slater values of a by com- 
puting the mean absolute error (MAE) in atomization 
energies of a set of 56 molecules (G2 set). We then de- 
termine the optimal value of a by minimizing the MAE 
for the G2 set of molecules. The calculations are per- 
formed for various basis sets in order to study the basis 
set dependence of the optimal a value. Our calculations 
show that the analytic HFS model with the optimal a 
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value performs better than the HF theory or the LDA 
and hence provides a computationally efficient scheme to 
study large systems at modest accuracy. Furthermore, 
by minimizing the MAE between the HF and the HFS 
total energies for atoms H through Kr, we find that best 
performance of the exchange functional in Eq.Q is ob- 
tained for a = 0.7091. 

Our calculations in the Slater-Roothaan (SR) method 
require using the Gaussian basis sets to fit the or- 
bitals and the Kohn-Sham potential. We have used the 
valence triple-C (TZ) 6-311G**^2I and the DGaussSS 
valence double-^— basis set (DZVP) for the orbitals. 
The s-type fitting bases are obtained by scaling the s- 
part of the orbital basisi^i For the non-zero angular 
momentum components the resolution-of-the-identity-J 
(RI-J) 25 and A22£ basis sets are used for the Kohn- 
Sham potential fitting. Thus, four sets 6311G**/RI-J, 
6311G*7A2, DZVP/RI-J, and DZVP/A2 of bases were 
used for optimizing the a value. The molecules were 
optimized using the Broyden-Fletcher-Goldfarb-Shanno 
(BFGS) algorithm^ The a minimization was performed 
using powerful Perl scripts that drive the analytic DFT 
code. 

Our first attempt to determine the optimal a is based 
on the atomic calculation. These calculations are nu- 
merical and therefore are free from the basis sets effects. 
Here we minimize the MAE in the HF and the HFS total 
energies for atoms H through Ar. The minimum occurs 
for a = 0.7267 with the MAE of 0.101 a.u. The opti- 
mal a value decreases slightly to 0.7091 when the target 
set is extended to include the atoms up to krypton. At 
this value the MAE is 0.33 au. These errors are an or- 
der of magnitude smaller than the MAE (2.38 a.u.) for 
the a — 2/3. The exchange functional in Eq. (Q with 
a = 0.7091 is therefore better approximation than the 
GKS functional, at least for the atomic systems. 

We now examine the performance of the HFS model 
for the GKS and Slater a values. The performance of the 
analytic HFS model in prediction of the atomization en- 
ergies for the G2 set of molecules is given in Tabled The 
computation of atomization energies is a stringent test 
for computational models and has been routinely used in 
the appraisal of the computational models. The G2 set 
of molecules used in the performance analysis is the set 
of 54 molecules in 56 electronic states due to Pople and 
coworkers^ and is often used for performance tests' It 
is apparent from the Table ^ that the errors are consid- 
erably smaller for the GKS value of a than for Slater a. 
The atomization energies are overestimated for Slater a 
with a mean error of 52.4 kcal/mol while for the GKS a 
the molecules are by and large underbound. The mini- 
mization of the MAE leads to the a value close to 0.7. 
The atomization energies for the G2 set at the optimal 
values are also shown in Tabled Going from the GKS a 
to the optimal a , the MAE reduces by about 2 kcal/mol 
while the mean errors decreases by 6-7 kcal/mol. The 
reduction in the mean error mainly occurs because it 
changes sign for more molecules. For the GKS a, the 



error is maximum (-61.49 kcal/mol) for the Si2Hg. It is 
also maximal for Si2Hg at the optimal a . CO2 is another 
molecule for which the error is comparable to this error 
at the optimal a . 

In order to investigate the role of basis sets on the 
optimal a value, we optimized the a for four different 
basis sets. The optimal a values are 0.70650, 0.69937, 
0.7032, and 0.698 for the 6311G**/RI-J, 6311G**/A2, 
DZVP/RI-J, and DZVP/A2, respectively. The MAE (in 
units of kcal/mol) at these optimal values are 13.5, 13.4, 
12.8, and 12.8, respectively. All optimal values are close 
to 0.7 and are effectively insensitive to the orbital basis 
sets. A small dependence on the fitting basis set is how- 
ever noticeable. The best performance is obtained for the 
DZVP/A2 basis set for which the MAE is 12.8 kcal/mol. 

We also carried out the performance test of the model 
in predicting the bond distances. For this purpose we 
selected all (15 in total) diatomic molecules belonging to 
the G2 set. Our results show that the MAE in bond 
distances also is smaller at the optimal a value. For 
the 6311G**/RI-J the MAE at optimal value is 0.019 A, 
0.013 A smaller than at the GKS a (0.032 A). The basis 
set effects show that the larger 6311G**/RI-J performs 
better than the DZVP/A2 basis. The MAE at their op- 
timal a values for these two bases are 0.019 and 0.048 
A, respectively. These are comparable or better than the 
LDA ( 0.024 A) or HF (0.028 A) errors^ 

In Table [H] we have summarized the results of present 
calculations. In comparison with the LDA or HF the an- 
alytic HFS model performance is significantly better. Its 
performance is even better than the SR-HF or the SR- 
Exact-Atomic models. These models are similar to the 
present one but make use of atom dependent a, which in 
case of SR-HF and SR-Exact- Atomic uses a values that 
give the HF atomic and the exact energies for atomic sys- 
tems, respectively. The overall improvement in the per- 
formance obtained here by minimizing the MAE in atom- 
ization energies also suggest that the SR model can also 
be similarly improved by multidimensional minimization 
of MAE in the a space. There are several density func- 
tional computational schemes that use the generalized 
gradient approximation (GGA), the hybrid GGA or meta 
GGA(See for example, Ref. ll9J). The accuracies of these 
models for the G2 set range from 3-8 kcal/mol, but to 
date they require numerical treatment. Although its an- 
alytic implementation is computationally most efficient, 
the optimal values can also be used in any existing den- 
sity functional code, albeit with some reduction in com- 
putational performance. It should also be borne in mind 
that the G2 set used in obtaining optimal a value con- 
tains small molecules consisting of atoms belonging to 
the first and second rows of the periodic table. 

We have also examined the performance of the ana- 
lytic HFS model for the extended G2 set containing 148 
molecules. Reoptimizing the a in order to minimize the 
MAE for this larger G2 set moves optimal a significantly 
far in the direction of the GKS's a value. The analysis 
of errors for individual molecules in this dataset shows 
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TABLE I: The atomization energies Do (kcal/mol) for the 56 set of molecules for two different basis sets. The two basis 
sets chosen are I: 6311G**/RIJ, H: DZVP/A2. The a values are Gaspar-Kohn-Sham's alpha (=0.66666667), Slater's alpha 
(=1.0000), and the optimal a for which mean absolute error is minimum. The last column contains the exact values. 





Basis I 


Basis I 


Basis I 


Basis II 


Basis II 


Basis II 




a 


0.666667 


1.00000 


0.70650 


0.666667 


1.00000 


0.69800 


Exact 


H 2 


81.56 


91.53 


82.76 


84.16 


94.13 


84.90 


103.50 


LiH 


33.44 


48.59 


34.86 


30.89 


41.90 


31.49 


56.00 


BeH 


46.62 


79.40 


50.09 


25.52 


46.23 


27.11 


46.90 


CH 


63.16 


74.57 


64.46 


64.95 


75.70 


65.78 


79.90 


CH 2 ( 3 Bi ) 


172.67 


242.62 


180.69 


175.69 


244.94 


181.41 


179.60 


CH 2 CM ) 


144.13 


187.45 


149.03 


148.50 


190.34 


151.92 


170.60 


CH 3 


269.96 


359.08 


280.32 


274.91 


362.50 


282.28 


289.20 


CH 4 


369.54 


484.63 


382.87 


375.72 


490.08 


385.38 


392.50 


NH 


62.39 


77.94 


64.13 


65.12 


79.59 


66.21 


79.00 


NH 2 


144.46 


186.28 


149.20 


150.51 


189.84 


153.69 


170.00 


NH 3 


247.52 


327.10 


256.68 


256.44 


332.45 


262.79 


276.70 


OH 


91.11 


116.00 


94.03 


94.17 


116.74 


96.03 


101.30 


H 2 


205.36 


281.07 


214.13 


213.20 


283.90 


219.17 


219.30 


HF 


129.94 


181.58 


136.01 


136.08 


183.42 


140.19 


135.20 


Li 2 


6.75 


7.78 


6.65 


5.74 


9.69 


5.66 


24.00 


LiF 


127.73 


200.55 


135.59 


118.91 


192.70 


124.75 


137.76 


C 2 H 2 


380.64 


508.89 


395.62 


375.39 


495.61 


385.98 


388.90 


C 2 H 4 


514.63 


688.39 


534.80 


519.43 


688.06 


534.03 


531.90 


C 2 He 


640.65 


858.65 


665.79 


649.78 


864.88 


668.10 


666.30 


CN 


175.85 


226.64 


180.87 


168.49 


210.50 


171.89 


176.60 


HCN 


291.80 


371.10 


300.99 


284.76 


354.24 


291.08 


301.80 


CO 


262.44 


327.45 


269.96 


251.66 


307.16 


257.01 


256.20 


HCO 


279.07 


370.60 


289.35 


274.94 


359.47 


282.37 


270.30 


H 2 CO 


357.03 


476.46 


370.63 


356.06 


469.07 


365.93 


357.20 


CH 3 OH 


468.42 


640.15 


487.97 


476.56 


643.49 


490.71 


480.80 


N 2 


206.16 


233.59 


209.47 


196.84 


212.00 


198.89 


225.10 


N 2 H 4 


371.99 


512.35 


387.70 


386.51 


519.20 


397.56 


405.40 


NO 


153.63 


186.81 


157.31 


146.96 


171.38 


149.52 


150.10 


o 2 


144.97 


191.03 


149.76 


141.18 


183.86 


145.09 


118.00 


H 2 2 


255.11 


356.29 


266.06 


262.57 


361.22 


270.64 


252.30 


F 2 


60.99 


86.35 


63.10 


60.70 


91.54 


63.08 


36.90 


C0 2 


413.28 


554.85 


429.25 


398.39 


524.83 


409.85 


381.90 


SiH 2 ( 1 A 1 ) 


116.69 


140.69 


119.57 


119.85 


146.07 


121.96 


144.40 


SiH 2 ( 3 Bx ) 


110.16 


149.90 


114.64 


112.65 


153.60 


115.87 


123.40 


SiH 3 


179.61 


229.59 


185.47 


183.25 


234.73 


187.39 


214.00 


SiH 4 


260.34 


326.68 


268.22 


264.75 


333.84 


270.39 


302.80 


PH 2 


119.12 


139.95 


121.62 


123.02 


147.19 


124.89 


144.70 


PH 3 


191.36 


229.18 


195.93 


197.19 


240.50 


200.61 


227.40 


H 2 S 


154.22 


188.60 


158.33 


158.87 


199.34 


162.23 


173.20 


HC1 


95.69 


118.24 


98.50 


96.37 


123.43 


98.65 


102.20 


Na 2 


5.47 


3.85 


5.20 


5.59 


4.66 


5.42 


16.60 


Si 2 


68.36 


90.14 


70.80 


68.42 


92.93 


70.48 


74.00 


P 2 


91.33 


101.21 


92.69 


91.67 


103.08 


92.73 


116.10 


S 2 


102.03 


137.80 


106.19 


103.03 


141.76 


106.46 


100.70 


Cl 2 


57.18 


92.08 


60.88 


58.37 


95.81 


61.29 


57.20 


NaCl 


80.52 


121.13 


84.89 


81.00 


127.66 


84.84 


97.50 


SiO 


182.85 


239.02 


189.13 


182.19 


239.66 


187.07 


190.50 


CS 


165.02 


202.06 


169.50 


163.97 


202.37 


167.64 


169.50 




1 98 41 


1 77 


IOO. UO 




1 83 3^ 
100.00 


1 37 30 

10 1 .ou 


1 93 


CIO 


69.56 


104.25 


72.98 


74.70 


111.20 


77.59 


63.30 


C1F 


69.97 


109.07 


73.98 


74.16 


116.73 


77.61 


60.30 


Si 2 H 6 


438.03 


562.19 


452.72 


445.26 


576.48 


456.04 


500.10 


CH3CI 


361.67 


484.96 


375.92 


366.70 


492.51 


377.46 


371.00 


H 3 CSH 


423.74 


560.88 


439.59 


432.27 


574.77 


444.35 


445.10 


HOC1 


156.11 


222.91 


163.21 


163.22 


234.11 


168.94 


156.30 


so 2 


250.06 


356.02 


261.65 


254.36 


360.56 


263.42 


254.00 


mean absolute 


15.9 


55.2 


13.5 


14.5 


55.3 


12.8 
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TABLE II: The mean absolute error (MAE) (kcal/mol) in 
the atomization energy of 56 molecules belonging to the G2 
set is compared within different models. The numbers for 
the SR-HF and SR- Exact-Atomic are for the Slater-Roothaan 
model with Hartree-Fock a values and the a values that give 
the exact-atomic (See text for more details). The results of 
the more complex PBE GGA functional are also included for 
comparison. 



Model 



Basis 



MAE 



Hartree-Fock theory 

LDA 

PBE 

SR-HF 

SR-HF 

SR- Exact-Atomic 
SR- Exact- Atomic 
HFS (Uniform a ) 
HFS (Uniform a ) 



6-311G**/RI- 

DZVP/A2 

6-311G**/RI- 

DZVP/A2 

6-311G**/RI- 

DZVP/A2 



78 Ref.28 

36 ReOO 
8 ReOO 

16 Ref.8 

16 Ref.8 

19 Ref£ 

18 Ref£ 

14 (Present work) 

13 (Present work) 



that this occurs due to the presence of a large percent of 
molecules containing fluorine in the extended G2 dataset. 
The errors for these molecules are lowered by decreasing 
the a value below 0.7. This is consistent with our ear- 
lier finding that the exact atomization of fluorine dimcr 
is obtained for much smaller a value of 0.3£ This again 
brings out the limitation of the uniform a HFS method 
and shows that the analytic SR method has a scope for 
improvement. It appears from the minimization of errors 
of the G2 and the extended G2 data sets and error anal- 
ysis, as well as from the minimization of the total total 



atomic energies that overall the value close to 0.7 is prob- 
ably the right choice for the optimal a in the uniform a 
calculations. 

To summarize, the performance appraisal of the ana- 
lytic Hartree-Fock-Slater method is carried out for vari- 
ous a values using the G2 database of 56 molecules. The 
a value that gives the best performance is determined 
by minimizing the mean absolute errors in the atomiza- 
tion energies of the G2 set of molecules. It is shown 
that the analytic HFS model performs better than the 
LDA or HF as well as the SR method that uses atom 
dependent a which give the exact HF or experimental 
atomic energies. Further, by minimizing the MAE in the 
HF and the HFS total energies it is shown that the lo- 
cal exchange functional performs significantly better for 
a = 0.7091 than the Gaspar-Kohn-Sham exchange func- 
tional. The MAE in former is an order of magnitude 
smaller than the MAE for the GKS exchange functional. 
The use of this exchange functional in more sophisticated 
GGAs could boost their performance considerably, and 
performance gain is already observed in case of Becke's 
exchange functional. 31 

Analytic DFT, even at this stage of development, is 
remarkably accurate. 
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* Electronic address: rzope@alchemy.nrl.navy.mil 
' Electronic address: dunlap@nrl.navy.mil 

1 K. S. Werpetinski and M. Cook, Phys. Rev. A 52, 3397 
(1995); J. Chem. Phys. 106, 7124 (1997). 

2 B. G. Johnson and M. J. Frisch, Chem. Phys. Lett. 216, 
133 (1993). 

3 J. A. Pople, Rev. Mod. Phys. 71, 1267 (1998). 

4 W. Kohn, Rev. Mod. Phys. 71, 1253 (1999). 

5 J. C. Slater, Phys. Rev. 81, 385 (1951). 

6 B. I. Dunlap, J. Phys. Chem. 107 10082 (2003). 

7 R. R. Zope and B. I. Dunlap, Chem. Phys. Lett 386, 403 
(2004). 

8 R. R. Zope and B. I. Dunlap (unpublished). 

9 P. A. M. Dirac, Proc. Cambridge Philos. Soc. 26, 376 
(1930). 

10 Caspar, R. Acta Phys. Hung. 3, 263 (1954). 

11 W. Kohn and L. J. Sham, Phys. Rev. 140 A1133 (1965). 

12 K. Schwarz, Phys. Rev. B 5, 2466 (1972). 

13 V. H. Smith and J. R. Sabin, J. Phys. B 11, 385 (1978). 

14 Connolly, J. W. D. In Modern Theoretical Chemistry; Se- 
gal, G. A., Ed.; Plenum: New York, 1977; Vol. 7, p 105. 

15 S. Fliszar S., N. Desmarais, M. Comeau, J. Mol. 
Struct. (THEOCHEM) 251, 83 (1991); Vauthier, E. C, 
Cosse-Barbi, A.; Blain, M.; Fliszar, S. J. Mol. Struct. 
(THEOCHEM) 492, 113 (1999). 

16 E. J. Baerends and P. Ros, Chem. Phys. 2, 52 (1973). 



17 E. J. Barends and P. Ros, Int. J. Quan. Chem. Symp. 12, 
169 (1978). 

18 O. Gunnarson, J, Harris, and R. O. Jones, Phys. Rev. B 
15, 3027 (1977). 

19 W. Koch and M. C. Holthausen, A Chemist's Guide to 
Density Functional Theory (Wiley- VCH: Weinheim, Ger- 
many, 2001) p 74. 

20 R. Krishnan et al. , J. Chem. Phys. 72, 650, (1980). 

21 A. D. McLean, G. S. J. Chandler, J. Chem. Phys. 72, 5639 
(1980). 

22 J. Andzelm, E. Wimmer, J. Phys. B, 172, 307 (1991); J. 
Chem. Phys. 96, 1280 (1992). 

23 N. Godbout et al. , Can. J. Chem. 70, 560 (1992). 

24 B. I. Dunlap, J. W. D. Connolly, and J. R. Sabin, J. Chem. 
Phys. 71, 3396; 4993 (1979). 

25 K. Eichkorn et al. , Theor. Chem. Acc. 97 (1997) 119. 

26 W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. 
Vetterling, Numerical Recipes The Art of Scientific Com- 
puting (Cambridge University Press: Cambridge, England, 
1986) p 309. 

27 J. A. Pople et al. , J. Chem. Phys. 90 5662 (1989); L. A. 
Curtiss et al. , J. Chem. Phys. 94 7221 (1991). 

28 A. D. Becke, J. Chem. Phys. 98 1372 (1993). 

29 B. G. Johnson, P. M. W. Gill, and J. A. Pople, J. Chem. 
Phys. 98 5612 (1993). 

30 M. Ernzerhof and G. E. Scuseria, J. Chem. Phys. 110 5029 



(1999). 

N. C. Handy and A. J. Cohen, Mol. Phys. 99, 403 (2001). 



